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ON-LINE COMBUSTIONLESS MEASUREMENT OF GASEOUS 
FUELS FED TO GAS CONSUMPTION DEVICES 

Background of the Invention 
Technical Field 

The present invention relates to methods for the 
combustionless measurement of the quality of gaseous 
fuels fed to gas consumption devices, and particularly 
natural gas consumption devices, as well as apparatus 
5 for the carrying out of the different variations of the 
process . 

Description of the Prior Art 

The heating value of a substance is of significant 

10 interest because it forms one basis for determining the 
commercial value of that substance as a fuel. Methods 
for measuring the quality of gaseous fuels to ascertain 
the amount of heat available therefrom are already being 
used in practice for numerous purposes. Recently, 

15 interest in and need for such measurements have 
increased considerably for various reasons. In 
industrial heating processes, it is frequently necessary 
to feed a well defined amount of heat per unit of time 
to a furnace in order to obtain optimum results. In 
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other cases it is desirable to optimize the consumption 
of fuel, i.e., to feed only the amount of heat actually 
required even if supplying a larger amount of heat does 
not adversely affect the process or product. For 
5 accounting purposes, billing on the basis of the amount 
of heat supplied has also been preferred to billing on a 

volume basis. 

An extensive need has arisen, therefore, for 
measurements of the quality of a gaseous fuel. 

10 Unfortunately, gas quality measurement is complicated by 
the fact that combustion gases, and particularly natural 
gases, are typically distributed together 
notwithstanding separate origin, composition and 
properties that differ to a greater or lesser extent 

15 from each other. Since processes and apparatus have 

proven themselves for use with such distributed gases of 
different composition, the on-line measurement of the 
quality of the gas or the quantity of heat available 
therefrom has gained increasing importance for the 

20 industrial use of gas and for accounting purposes. 

Known methods for evaluating gas quality are in 
most cases not readily adaptable to these uses due to 
technical reasons or because of cost considerations. 
For example, one conventional method for measuring 

25 heating value known in the art comprises combustion 
calorimetry. This process involves the burning of a 
partial stream of the combustible gas with an open flame 
or with a catalyst and measuring the heat produced. The 
necessity of burning a measured partial stream of the 

30 gas in order to determine its heating value, as known 
from experience, requires frequent maintenance of the 
apparatus, since a flame can change due to deposits of 
combustion residues or because a combustion catalyst 



WO 93/08457 



PCT/US92/09150 



gradually declines in effectiveness. The required 
accuracy of measurements which serve, for example, for 
billing purposes can only be obtained if these apparatus 
are operated under well-defined, controlled conditions, 
5 preferably in an air-conditioned chamber, which is 
obviously expensive* 

Other known non-combustion methods for continuously 
analyzing a stream of gas include gas chromatography and 
mass spectrometry. Gas chromatography and mass 

10 spectrometry are techniques for separating and 

identifying each constituent of the gas and measuring 
the relative concentration thereof. Knowing the heating 
value of each constituent of a mixture, the total 
heating value may then be computed. Unfortunately, 

15 these methods require a large expenditure of measurement 
and control devices to implement. 

The same techniques are also currently used to 
determine other parameters representative of the quality 
of a distributed fuel gas, such as density and percent 

20 concentration of inert gases therein. A principle use 
for gas density determination is in the operation of an 
orifice flow meter, while percent concentration of inert 
gases, for example, nitrogen, N 2 , carbon dioxide, C0 2 , 
and oxygen, 0 2 , is used to determine gas 

25 pumping/ transportat ion cost or for subsequent regulation 
of a combustion process. 

To summarize, most, if not all, presently known 
techniques for determining the quality of a fuel gas, 
such as heating value, density or percent concentration 

30 of inert gases, have one or more drawbacks associated 
therewith, including: requiring trained personnel to 
operate, producing time delayed results, lacking 
repeatability, destroying the sample, being cumbersome 
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or expensive to implement and lacking sufficient 
accuracy due to an inability to completely distinguish 
constituents. Therefore, there exists a genuine need in 
the art for a novel approach to the measurement of the 
5 quality of a fuel gas which is accurate, reliable and 
inexpensive in implementation. 

Summary of th e Invention 
The primary object of the present invention is to 
10 provide a process and apparatus capable of determining 
the quality of a combustible gas without the combustion 
of the gas and in a simpler manner than with previously 
known combustionless methods. 

Another object of the present invention is to 
15 provide such a process and apparatus which are more 
readily adaptable to field installation than 
conventional techniques for effecting such measurements. 

Yet another object of the present invention is to 
provide such a process and apparatus which are capable 
20 of determining the heating value, density or percent 
concentration of gas inerts of the fuel gas. 

A further object of the present invention is to 
provide such a process and apparatus which are capable 
of determining one or more of the heating value, density 
25 and percent concentration of inerts within the fuel gas 
from the same readily measured gas parameters. 

The foregoing and other objects are accomplished in 
accordance with the present invention in part through 
the discovery of an empirical formula correlating 
30 certain readily measurable gas parameters with the 
desired measurement, i.e. r - heat content, density or 
percent concentration of gas inerts. In their broadest 
aspect, the process and apparatus recited herein utilize 
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a formula of the form: 

M - a 0 + b^," 1 + bp** + ... + + c^ 02 + ... (1) 

+ d iSti Pl + + ••• + e Ai Ul + **tz* + ••• 

where: 

5 a 0 , b,, b 2 , c t , c^, d 1# d 2 , e 1# % « constants; 

nl, n2, ml, m2, pi, p2 f ul, u2 = exponents; 
n ti = viscosity at various temperatures, ti; 
= thermal conductivities at various 
temperatures, ti; 
10 c pt - = specific heat at various temperatures, ti; 

and 

A ti = optical absorption at various temperatures, 
ti. 

Those skilled in the art will appreciate that 
15 additional, readily measured gas parameters, such as 

speed of sound, may also be incorporated into the above 
formula, provided the overall novel form of the equation 
is maintained. Also, one or more measured terms, for 
example, viscosity, specific heat and/ or optical 
20 absorption, may be omitted from the formula if 

unnecessary to attaining a desired accuracy level. At 
least two different terms are believed necessary, 
however. 

By way of example, in one specific embodiment 
25 applicants empirical formula is expressed as: 

M = a c + b.ffn) 01 + c 1 f l (k t1 ,k t2 ) B1 + 0^(1^,**) 012 (2) 
where: 

M = heat content, density or percent gas inerts 
calculation; 
30 a c , b 1# c, & c 2 ■ constants; 

ol, ml & m2 » exponents; 
n = viscosity; 

= thermal conductivity at a first temperature; 
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and 

Jc^ = thermal conductivity at a second temperature. 
Prior to ascertaining a desired value (i.e., heat 
content, density or percent concentration of gas 
5 inerts), applicants method requires the steps of: 
conducting at least a partial stream of the fuel gas 
through a sensor chamber having a plurality of sensors 
in contact with a fuel gas; generating a first 
electrical signal at one of the plurality of sensors, 

10 the first electrical signal being representative of a 
first fuel gas quality, the first fuel gas quality 
comprising one of thermal conductivity, specific heat, 
viscosity and optical absorption; conducting the first 
electrical signal to a computing means; generating a 

15 second electrical signal at one of the plurality of 
sensors, the second electrical signal being 
representative of a second fuel gas quality, the second 
fuel gas quality comprising one of thermal conductivity, 
specific heat, viscosity and optical absorption, the 

20 second fuel gas quality comprising a different one of 
the fuel gas qualities than the first gas quality; 
conducting the second electrical signal to the computer 
means; and finally, using said computing means to derive 
a signal for at least one of measurement utilizing the 

25 first and second electrical signals as a measure for at 
least one of the heat content, density and percent gas 
inerts of the fuel gas according to formula: 

M = a 0 + b,*" 1 +.c iY p 1 

where: 

30 n = said at least one of said fuel gas heat 

content, density and percent concentration of 
gas inerts, 
a c , b v c r , - constants, 
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ml, pi, = exponents, 

x = signal representative of the first fuel gas 
quality, and 

y = signal representative of the second fuel gas 
5 quality. 

Specific values for the constants and exponents, which 
have been defined by applying linear progression 
analysis to experimental test results, are provided. 

In another aspect, the present invention comprises 

10 a corresponding apparatus for the combust ionl ess 
measurement of fuel gas. The apparatus includes a 
sensor chamber having a plurality of sensors therein and 
conducting means for moving at least a partial stream of 
the fuel gas through the sensor chamber such that the 

15 gas is in contact with the plurality of sensors. First 
generating means is provided for producing a first 
electrical signal at one of the plurality of sensors. 
The first electrical signal is representative of a first 
fuel gas quality, which comprises one of thermal 

20 conductivity, specific heat, viscosity and optical 

absorption of the fuel gas. Transferring means conducts 
the first electrical signal to a computer for 
processing. A second generating means is also included 
for generating a second electrical signal at one of the 

25 plurality of sensors. The second electrical signal is 
representative of a second fuel gas quality which is 
different from the first fuel gas quality. The second 
fuel gas quality comprises one of thermal conductivity, 
specific heat, viscosity and optical absorption of the 

30 fuel. Again, transferring means conducts the second 
electrical signal to the computer for processing. The 
computer is then used to periodically derive a signal 
for at least one of measurement and regulation using the 
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f irst and second generated electrical signals as a 
measure for at least one of the heat content, density 
and percent gas inerts of the fuel gas according to the 
formula: 

5 M - a 0 + b,^ 1 + c^ 1 

where: 

H - said one of said fuel gas heat content, 
density and percent concentration of gas 
inerts, 

10 a of b,, c t , « constants, 

ml, pi, ■ exponents, 

x = signal representative of the first fuel gas 
quality, and 

y = signal representative of the second fuel gas 

15 quality. 

In a preferred form, the plurality of sensors comprise 
microsensors which are arranged within the chambers such 
that substantially zero gas flow is encountered thereby. 
Additional embodiments of the process and apparatus of 

20 the present invention are described below. 

Brief Description of the Drawings 
The subject matter which is regarded as the 
invention is particularly pointed out and distinctly 
25 claimed in the concluding portion of the specification. 
The invention, however, both as to organization and 
method of practice, together with further objects and 
advantages thereof, may best be understood by reference 
to the following description taken in connection with 
30 the accompanying drawings in which: 

Fig. 1 is an operational overview of one process 
embodiment of the present invention; 

Fig. 2 is a schematic illustration of one 
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embodiment of the apparatus of the present invention; 

Fig. 3 is a block diagram of the analog circuit 
board depicted in Fig. 1; 

Figs. 4(a), 4(b), and 4(c), represent several 
5 heater/sensor configurations of microbridge systems in 
accordance with applicants preferred implementation of 
the present invention; 

Fig. 5 is a schematic representation of sensor 
time/temperature response curves according to a heater 
10 pulse; 

Fig. 6 is a partial schematic and block diagram of 
a circuit for use with a sensor as depicted in Fig. 4(b) 
in accordance with the preferred implementation; 

Fig. 7 is a more detailed circuit schematic with 
15 reference to Fig. 4(c); and 

Fig. 8 is a schematic block diagram of the 
preferred thermal conductivity sensor embodiment 
including calibration and use functions. 

20 Detailed Description of the Invention 

As briefly set forth above, central to the manner 
of technical action which forms the subject matter of 
this invention, as defined in the attached process and 
apparatus claims, is the surprising discovery iAat heat 

25 content, density and percent concentration of gas inerts 
may each be readily and accurately determined from an 
empirical expression, for example, of the form: 

,M = a 0 + b 1 f(n)° 1 + c l f 1 (k t1 ,Jc t2 ) B1 + c^f 2(^,1^) 02 (2) 
30 where: 

M = heat content, density or percent gas inerts 

calculation; 
a Q , b 1# c, & c 2 = constants; 
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10 



ol, ml & m2 = exponents; 
n = gas viscosity; 

k^ = gas thermal conductivity at a first 

temperature, tl; and 
Jc^ = gas thermal conductivity at a second 35 
temperature, t2. 
(Units for the various measured parameters and 
calculations are set forth in Table II below.) Applicant 
has also discovered that several other readily measured 
parameters characteristic of the fuel gas, such as 
specific heat, c p , and optical absorption, A, may be 
determined and used to supplement, or substitute for, 
the viscosity, n, and thermal conductivity, k, variables 
of equation (2) (e.g., see equation (1) above and 
15 equation (8) below). Notwithstanding the utilization of 
different gas parameter combinations, however, the 
overall form of equations (1), (2), (8), etc., is 
maintained. (As set forth, equation (2) in part 
expresses gas thermal conductivity as a function of two 
20 readings k,, and k^ (e.g., see equations (3) - (7) 
below) . One skilled in the mathematics art could 
separate these terms if desired to attain a formula more 
clearly of the form of equation (1).) All the gas 
parameters discussed herein are readily measurable using 
25 existing technology, but preferred implementing 

processes and apparatus are described, and claimed, 
below. 

Specific algorithm examples are provided herein for 
determining the quality of natural gases. However, 
30 applicant believes that the empirical formula is equally 
applicable to other types of fueJ^ gases and that one of 
ordinary skill can derive the necessary specific 
equations therefor from the information provided herein. 
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Figure 1 depicts one operational overview of a 
method of the present invention implemented using the 
algorithm of equation (2) . Initially, a portion of the 
fuel gas must be diverted through a sensor chamber 
5 (discussed below) , 10 "Divert Portion Of Flowing Fuel 
Gas Through Sensor Chamber." First and second 
measurements 12 "Measure Fuel Gas Viscosity At A First 
Sensor" and 14 "At A First Temp. , Measure Fuel Gas 
Thermal Conductivity At A Second Sensor," respectively , 

10 are taken. Thereafter the temperature of the fuel gas 
at the second sensor is increased, 16 "Increase Temp. 
Of Fuel Gas Adjacent Second Sensor," and the fuel gas 
thermal conductivity at the increased temperature is 
determined, 18 "Measure At Increased Temp. Fuel Gas 

15 Conductivity At Second Sensor." In the preferred 

embodiment, an increase in fuel gas temperature adjacent 
the second sensor is readily attained by increasing the 
current though a microbridge structure (described below) 
used to obtain the thermal conductivity measurements. 

20 The two temperatures at which gas thermal conductivity 

is measured should be selected to optimize the signal to 
noise ratio of the resultant microbridge output. If too 
low a temperature is initially selected as tl, the 
difference signal is not strong enough, and if too high 

25 a second temperature is chosen for t2, the operational 
lifetime of the sensor will be short. Approximately 
70° C and 120° C are believed to be examples of 
acceptable first temperature, tl, and second 
temperature, t2, values for ascertaining gas thermal 

30 conductivity with the arrangement described herein. 

Next, the method requires correction of measured 
viscosity and thermal conductivity values to account for 
influences of gas temperature, gas pressure and circuit 
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temperature changes, 20 "Correct Fuel Gas Viscosity and 
Thermal Conductivity Measurements Based Upon Fuel Gas 
Temp, and Pressure, and Sensor Circuit Temp." Viscosity 
and thermal conductivity values can also be converted to 
standardized readings using the measured absolute 
pressure and temperature of the fuel gas. Equation (2) 
is then used to determine the desired heating value, 
density or percent concentration of gas inerts value for 
the fuel gas, 22 "Calculate Fuel Gas Heating Value, 
Density and/or Percent Concentration of Inert Gases." 
This value is then either stored or displayed, 24 "Store 
or Display Fuel Gas Quality Calculation," and thereafter 
flow returns to step 12, i.e., subsequent a time delay, 
26 "Time Delay." Time delay 26 must be sufficient for 
15 the thermal conductivity sensor to return to first 
temperature level tl, i.e., the temperature of the 
sensor chamber. With a microbridge sensor configuration 
as described below, this only requires a few 
milliseconds. 

20 By way of more specific example, formula (2) is 

first used in the process and apparatus of the present 
invention to combustionlessly determine the heating 
value or heat content of a fuel gas, such as a natural 
or synthetic fuel. 
25 With the use of any commercially available linear 

progression analysis program, one specific equation for 
deriving the heat content of a natural gas, which takes 
the form of equation (2), comprises: 

H e = 3643.53 + 1050.71(102/n) 3 (3) 
30 - 7.602211^ - 2294.2 (k^/k^) 

where: 

H c = gas heat content; 
n = viscosity; 
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k t1 = gas thermal conductivity at a first 

temperature, tl; and 
k^ = gas thermal conductivity at a second 
temperature, t2. 
5 Although preferred techniques for measuring viscosity 
and thermal conductivity are described below, both 
parameters can be readily determined by one of ordinary 
skill in the art using presently available technology. 
For example, viscosity is measurable by determining the 

10 pressure drop across a capillary through which a known 
volume of gas is pumped with a positive displacement 
pump and thermal conductivity can be determined using a 
conventional wheatstone bridge circuit. From extensive 
testing, applicant has determined that utilization of 

15 equation (3) to determine heat content of a natural gas 
results in a maximum error of .06696 MJ/m 3 , with a 
standard error of .01831 MJ/m 3 . This result is 
considered well within acceptable error limits imposed 
by the industry. 

20 Once obtained, derived heat content values are 

indicated or transmitted to recording instruments or 
given off as control pulses depending upon the 
measurement and/or regulation information required for a 
particular application. The process and arrangement are 

25 also suitable for installation in measurement stations 
of high-pressure, long distance gas transmission (pipe) 
lines in which the (heat) quantity of flow is 
continuously recorded. 

An important benefit flowing from utilization of a 

30 formula of the type of equations (1) and (2) is the 
synergistic benefits obtained therefrom when 
incorporated into a process and arrangement as claimed 
herein. Namely, any one of heat content, density and 
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percent concentration of gas inerts of the fuel gas may 
be readily computed using the same basic measurements, 
e.g., viscosity and thermal conductivity, and the same 
basic algorithm type, e.g., equation (1) , (2), or C8) • 
5 This is because each of these values comprises a measure 
of the "quality" of the fuel gas. 

Fuel gas density determination is important to 
computing the orifice coefficient of the gas, which is 
necessary for operation of a typical orifice flow meter, 

10 along with determining gas compressibility factors and 
densitometer calculations. By using formula (2) , the 
present invention provides an inexpensive and accurate 
means for ascertaining density of a natural gas. 
Applying linear progression analysis to equation (2) and 

15 independently derived density readings, one specific 
formula for determining the density of natural gas is 
obtained: 

p = 4.3077 + . 22937 (102/n) 3 (4) 
- .012094k t1 - 2.2881 (kjj/JCflj 

20 where: 

p = gas density; 
n - gas viscosity? 

= gas thermal conductivity at a first 
temperature, tl; and 
25 = gas thermal conductivity at a second 

temperature, t2. 
As noted, quantifying of percent concentration of 
inert gas has traditionally been accomplished with 
sophisticated, costly and labor intensive means, such as 
30 gas chromotagraphy equipment. Accurate determination of 
inert gases like nitrogen, N 2 , carbon dioxide, C0 2 , and 
oxygen, 0 2 in fuel gas, and particularly natural gas, is 
important because these gases reduce the heating value 
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of the fuel gas, cause the pumping/ transportation cost 
per unit of gas energy to be increased, and reduce the 
supplier's revenue per delivered unit volume of gas, 
etc. Through extensive experimentation and testing, 
5 equations (5) -(7) below have been identified as being 
preferred specific forms of formula (2) for various 
combinations of inert gases within natural gas. Again, 
the constants and exponents were determined by applying 
linear progression analysis to independently measured 
10 results. 

Xii^+cqp " 288.69 - 23.818/n 3 (5) 
- .595751^, - 173.65(^2/^) 

15 XII2402 = 464.65 + 9.8185/n 3 - .42180^ (6) 

- 356.18(k t2 /k tl ) 

X-^ = -175.96 - 33.636/n 3 - .173951^ (7) 
+ 182.52 (ktf/ktt) 

20 where: 

X N2*o2*co2 - % concentration of N 2 , 0 2 and CO z in fuel 
gas, 

*N2+o2 = % concentration of N 2 and 0 2 in fuel gas, 
= % concentration of C0 2 in fuel gas, 
25 n = fuel gas viscosity, 

k t1 « gas thermal conductivity at a first 

temperature, t, and 
k t2 = gas thermal conductivity at a second 
temperature, t^. 
30 As noted briefly above, specific heat and optical 

absorption comprise two additional readily measurable 
fuel gas characteristics which may be added to formula 
(2), e.g., to improve accuracy of the resultant 
computations, or to substitute for one or the other of 
35 the viscosity and thermal conductivity variables. A 
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generalized egression of the discovered algorithm is 
equation (1) above. Table I lists several specific 
algorithms for calculating heating value of natural gas 
which have been derived from formula (1) and their 
5 measured accuracy levels. The formulas of Table I are 
considered merely illustrative of various parameter 
combinations possible pursuant to the invention 
described herein. From the present description, other 
specific formulas for determining heating value, 
10 density and/or percent concentration inerts can be 
readily derived by a person skilled in the art. The 
appended process and arrangement claims are operable 
with any such formula of a form derived from equation 
(l). 

15 m addition to the above, applicant has discovered 

that another particularly preferred expression for 
determining heat content of a natural gas comprises: 

H = -1287.7 + 808.700C • 73846 . ^./ 8 > 

20 ° - 1,048, 800k* f ' 7<2 -.00090189 (Mn) 1 ' 7514 

where: 

c p = specific heat of the fuel gas; 
k = gas thermal conductivity; 

Mh = (molecular weight of thegas) • (viscosity of 

25 the gas). 

The term molecular weight times viscosity, Mn, or its 
alternate expression density (p) times viscosity, pn, is 
capable of being determined using a combination of 
available technologies. However, this quantity is 

30 preferably ascertained with the novel process and 
apparatus described in a copending U.S. patent 
application entitled "Multiple Gas^.Broperty Sensor", 
Serial No. 07/781,770 (see below for further 
discussion) . 
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The implementing apparatus of the invention will 
now be described in greater detail with reference to 
accompanying Figures 2-8. 

Referring first to Fig. 2, a structural 
implementation of one embodiment of the invention is 
shown. A fuel gas, such as natural gas, flows in main 
gas line 30 in the direction of arrows 32. A secondary 
gas inlet pipe 34 is provided for diverting a portion of 
the gas flow to a sensor chamber 36. The diverted gas 
is returned from chamber 36 to main line 30 via a gas 
outlet pipe 38, which connects to line 30 downstream 
from sensor chamber 36. A constriction 40 within main 
line 30 creates a pressure differential which forces a 
portion of the gas flow through sensor chamber 36. 
15 Radial dimensions for line 30 and inlet and outlet pipes 
34 and 38 may vary, but as an example, main line 30 is 
typically one to three inches in diameter, while pipes 
34 & 38 are approximately one-quarter inch in diameter. 
As shown, chamber 36 is located above main line 30, 
20 which minimizes condensation effects within the chamber. 
Sensor chamber 36 preferably comprises a metallic block 
into which inlet pipe 34 passes and from which outlet 
pipe 38 extends. Within the block is an airflow sensor 
42, a first pressure sensor 44 and a plurality of 
25 sensors 46, 48, 50 and 52, which, in the embodiment 

illustrated, are separated from gas flow by several flow 
blocking screens 54 such that substantially zero flow 
effects are encountered within area 56 of chamber 36 
adjacent said sensors. As discussed below, sensors 46, 
48, 50 and 52 preferably comprise microsensors because 
of their relative inexpensive cost and high accuracy. 
However, since these sensors are flow sensitive, screens 
54 are needed to create area 56 of minimal flow. 



30 
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Screens 54, e.g., manufactured of a fine wire mesh, are 
designed to prevent microscopic gas flow to area 56. A 
plurality of screens is used to ensure that fuel gas 
essentially only diffuses to area 56. A certain amount 
5 of turbulence exists within chamber 36 below screens 54. 
Sensors 44, 46, 48, 50 & 52 are electrically connected 
to an analog circuit board 58 via lead lines 60, which 
are electrically insulated from one another to prevent 
shorting. Leads 60 serve to connect the inside of 
10 sensor chamber 36 to the outside of the chamber. 

Preferably, insulation (not shown) surrounds chamber 36 
to facilitate maintenance of a constant temperature 
therein. 

In one embodiment, sensor 46 comprises a pressure 

15 sensor; sensor 48, a gas temperature sensor; sensor 50, 
a thermal conductivity sensor; and sensor 52, a 
viscosity sensor. In general, sensors 50 & 52 provide 
the information needed for calculation of heating value, 
density or percent gas inerts, e.g., using formula (2), 

20 while sensors 44, 46 & 48 allow for correction of sensed 
thermal conductivity and viscosity values for influences 
of pressure and temperature in a well known manner. In 
an alternate embodiment described below, only one 
sensor, e.g. sensor 48, is required to obtain the 

25 parameters necessary (i.e., thermal conductivity and 
specific heat) to determine heating value, density 
and/or percent gas inerts. 

Flow sensor circuitry 62 is provided to translate 
flow sensed at 42 into digitally readable signals. 

30 Sensor 42 and associated circuitry 62 are considered 
optional to implementation of the present invention • 
Digital processor 64 comprises any commercially 
available central processing unit or microcomputer. 
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Once the appropriate heating value, density and/or 
percent gas inerts determination has been realized, 
using one of the above formulas, the determination is 
outputted to a communications interface €6, which 
5 includes display means 68 and storage means 70 capable 
of presenting for viewing and/or storing determined 
values for subsequent use. If desired, a 
regulation/control signal is outputted via line 72. 
Analog circuit board 58 is depicted in greater 

10 detail in Fig. 3. 

As shown, circuit 58 includes pressure sensor 
circuitry 78 & 80, gas temperature sensor circuitry 82, 
circuit board temperature sensor circuitry 84 , thermal 
conductivity microbridge sensor circuitry 86, and its 

15 microbridge driver 90, and viscosity sensor circuitry 
88. Each sensor circuitry is electrically connected to 
the appropriate sensor and to an analog to digital 
converter 92. From converter 92, signals flow to 
processor 64 (Figure 2) . 

20 in the embodiment depicted, pressure sensor 

circuitry 78 is connected to sensor 44, and pressure 
sensor circuitry 80 is connected to sensor 46. 
Circuitry 78 & 80 and sensors 44 & 46 each comprise any 
commercially available pressure sensor, such as that 

25 manufactured by Honeywell, Inc. of Bloomington, 
Minnesota and marketed under Model No. ST3000. 
Temperature sensor circuitry 82 & 84 similarly comprise 
any well-known resistance temperature measurement 
apparatus such as a wheatstone bridge configuration 

30 wherein a change in the resistance branch exposed to the 
unknown temperature produces an unbalance proportional 
to said temperature. Circuitry 84 is connected to a 
sensor 85 which is positioned on analog circuit board 
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58. Circuitry 78, 80 , 82 & 84 and their associated 
sensors are necessary in order to account for sensor and 
circuit non-linearities and influence of electronic 
temperature changes, e.g., on resultant thermal 
5 conductivity and viscosity measurements. Also, actual 
measurements can be converted by one of ordinary skill 
in the art from measured gas temperature and pressure 
values to corresponding values at standard temperature 
and pressure. 

10 As noted above, there are numerous known techniques 

for determining thermal conductivity (and specific heat) 
of the fuel gas which could be implemented as sensor 
circuitry 86. However, applicant believes that a 
preferable approach is described in a recently issued 

15 patent entitled, "Measurement of Thermal Conductivity 
and Specific Heat," U.S. Patent No. 4,944,035, the 
entirety of which is hereby incorporated herein by 
reference. This approach will now be described in 
detail with reference to Figs. 4 (a) -8. 

20 With respect to measuring thermal conductivity in 

fluids, various typ6s of detectors have been used. This 
includes resistance bridge type sensors . _ One such 
device is described in U.S. Patent 4,735,082 in which 
thermal conductivity is detected using a wheatstone 

25 bridge technique in which a filament in one diagonal of 
the bridge is placed or positioned in a cavity through 
which the sample gas of interest is passed. The 
filament is used to introduce a series of amounts of 
thermal energy into the fluid of interest at alternating 

30 levels by varying the input voltage which, are, in turn, 
detected at the other diagonal as voltage difference 
signals. Integration of the changes of the value of the 
successive stream of signals yields a signal indicative 
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of the heat dissipation through the fluid, and thus, the 
thermal conductivity of the fluid. 

Further to the measurement of thermally induced 
changes in electrical resistance, as will be discussed 
in greater detail below, recently very small and very 
accurate -microbridge" semiconductor chip sensors have 
been described in which etched semiconductor 
"microbridges" are used as condition or flow sensors. 
Such sensors might include, for example, a pair of thin 
film sensors around a thin film heater. Semiconductor 
chip sensors of the class described are treated in a 
more detailed manner in one or more patents, such as 
U.S. Patent Nos. 4,478,076, 4,478,077, 4,501,144, 
4,651,564 and 4,683,159. 
15 It is apparent, however, that it has often been 

necessary to address the measurement of specific heat 
c , and thermal conductance, k, of a fluid of interest 
with separate and distinct devices. Not only is this 
quite expensive, it also has other drawbacks. For 
20 example, the necessity of separate instruments to 

determine specific heat and thermal conductivity may not 
allow the data consistency and accuracy needed for 
useful fluid process stream (gas or liquid) 
characterization. Further," the required degree of 
25 correlation may not be present. Because the 

determination of heat content as contemplated herein 
depends on both measurements, e.g., see equation (8) and 
Table I, this takes on even more importance. 

The referenced application overcomes the many 
30 disadvantages associated with the determination of both 
specific heat, c p , and thermal conductivity, k, by 
providing simple techniques which allow accurate 
determination of both properties in a sample of interest 
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using a single sensing system. The approach 
contemplates generating an energy or temperature pulse 
in one or more heater elements disposed in and closely 
coupled to the fluid medium (gas or liquid) of interest. 
5 Characteristic values of k and c p , of the fluid of 

interest then cause corresponding changes in the time 
variable temperature response of the heater to the 
pulse. Under relatively static sample flow conditions 
this, in turn, induces corresponding changes in the 

10 time-variable response of one or more temperature 

responsive sensors coupled to the heater principally via 
the fluid medium of interest. 

The thermal pulse of a source need be only of 
sufficient duration that the heater achieves a 

15 substantially steady-state tempera ture for a short time. 
This pulse produces both steady-state and transient 
conditions at the sensor. Thermal conductivity, k, and 
specific heat, c p , can be sensed within the same sensed 
thermal pulse by using the steady-state temperature 

20 plateau to determine k, which is then used with the rate 
of change of temperature in the transient condition to 
determine c p . Both values then provide input to the 
determination of heat content, density or percent 
concentration of gas inerts. 

25 The microbridge semiconductor chip sensor 

contemplated, for example, in certain embodiments 
preferred for the implementation of the 
thermoconductivity sensor may resemble the form of one 
or more of the microbridge systems illustrated in the 

30 patents identified above. Such a system is exemplified 
by the figures provided with Patent No. 4,501,144. 
Reference should be made to said patent for a better 
understanding of the discussion to follow. While the 
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present discussion is believed sufficient, to the extent 
necessary, additional material contained in the 
microbridge related patents cited is deemed to be 
incorporated herein by reference. 
5 Now with reference to the implementation of the 

thermal conductivity sensor, Figs. 4(a), 4(b) and 4(c), 
depict three slightly differing embodiments or 
configurations representative in terms of number and 
arrangement of the heaters and sensors which can be 

10 used in this structure. In Figure 4(a), all of the 

elements 122, 124 and 126 are used as heaters. Figure 
4(b) is an embodiment in which the thin film element 
126 acts as heater and elements 122 and 124 act as 
sensors. The embodiment of Figure 4(c), represents the 

15 preferred arrangement in which the element 122 acts as 
heater and element 124 acts as sensor. The effective 
gap and thus the thermal isolation between heater and 
sensor is desirably wider in the embodiment of Figure 
4(c). 

20 In the implementation of the applicant's preferred 

thermal conductivity sensor, particular attention is 
directed to (1) setting specific temperature markers in 
the sensor to determine the time periods needed for 
achieving the corresponding temperature changes, (2) 

25 using temperature sensors which are physically separated 
from the heater so that the direct influence of the 
heater and heat conducted to the sensor other than via 
the fluid of interest is reduced, and (3) using a pulse 
which reaches at least a momentary steady-state plateau 

30 to determine k, which then is used with the transient 
measure to determine c p . 

Figure 5 graphically depicts a square wave 
electrical energy pulse 130 to the heater 126 which 
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results in quasi square wave heat pulses released by the 
heater. These, in turn, result in reactive curves as at 
131, 132 and 133 at the sensor which vary as described 
below. The pulse applied to the heater, for example, 
5 may have a height of about 4 volts with a pulse width of 
100 ms. Since the heater is closely coupled through the 
fluid medium to the sensors, the family of curves 131, 
132 and 133 resembles the shape of the input pulse 130. 
They show the heat response in the sensors 122 and 124. 

10 The curves generally include beginning and ending 

transient portions flanking a relatively steady-state 
central portion. The quick response of the sensor 
allows a relatively long steady-state to exist even with 
a pulse of 100 ms. Of course, the curves are affected 

15 by factors such as pressure and temperature as they 
influence the effective thermal conductivity and 
specific heat of the particular fluid of interest. 

Heat flowing from the heater element or elements to 
the sensor element or elements is conducted both through 

20 the fluid and through the solid semiconductor element 
support substrate or the like. It is advantageous with 
respect to the measurement of k or c p of the fluid of 
interest that the amount of heat reaching the sensor 
through the solid connections be minimized so that 

25 substantially all the measured thermal effect is 
generated via the fluid of interest. 

With respect to the transfer of heat to the 
sensor (s) , some background information regarding the 
propagation of heat or temperature waves is presented. 

30 The speed of propagation, v, of a one dimensional wave 
(if it features an exponential decay profile) is 
constant and given by the expression: 

v = D T /ce = (D T /b) 0 - 5 , (9) 
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where: 

a = exponential decay constant, 
b m rise time constant at a fixed 
location, and 
5 D T = thermal diffusivity. 

A complete list of nomenclature and subscripts with 
units appears in Table II, below. D T is related to k 
and c by the expression: 

D T = k/c p < 10 > 
10 s D T , therefore, if known, may be a key to 

obtaining c p . The rise time constant, b, was measured 
to be about 4 msec. For typical gases, D T ranges from 
1.7 cm 2 /s for He to .054 cm 2 /s for CjH,,. Metals exhibit 
high values such as 1.7, 1.1 and .18 cm 2 /s respectively 
15 for Ag, Cu and Fe. Insulators, however, are even lower 
than the gases at .004 cm 2 /s for glass and .0068 cm 2 for 
Si 3 N 4 which, as discussed above, is a good insulator. 
The propagation speed, v, in a typical gas sample then 
is about (1/0. 004) 0 5 = 15 cm/s. This compares with 
20 (0. 0068/0. 004) 0 - 5 - 1.3 cm/s for Si 3 N 4 , assuming that the 
same rise time constant of about 4 is is applicable to 
both the one measured in the Si 3 N 4 and the actual one in 
the gas. 

The effect is that the influence of the temperature 
25 wave propagating from one thin film strip, that is, the 
heater, to a second thin film strip, the sensor, both 
being embedded in a membrane of Si 3 N 4 , since it reduces 
the contribution of heat flow through the solid media. 
This is beneficial to the accuracy of the system. 
30 Typical microbridge embodiments are illustrated by 

Figures 4(a) - 4(c) . They will now be explained in 

greater detail. 

The configuration of Figure 4(a) involves using the 



WO 93/08457 



PCT/US92/09150 



-27- 



same microresistance 122, 124, 126 for the heating pulse 
and the sensing task. In this embodiment, the resistive 
heater-sensor element may be one leg of a conventional 
resistive wheatstone bridge in a control circuit. 
5 TABLE II - NOMENCLATURE 

Symbol Units 
a Exponential Decay Constant cm 

a^ Constant 
A Optical Absorption 

10 B Area of Heat Transfer to Microbridge cm 2 

or to Gas 

b Rise Time Constant at a Fixed Location °C/s 

c p Specific Heat cal/(cm3°C) 

D t Thermal Diffusivity, D T = k/c p cm 2 /s 

15 H c Heat Content of Gas MJ/m 3 

k Thermal Conductivity cal/(sm°C) 

L Length of Thermal Conductance Path cm 

in Gas or Solid 
M Molecular Weight of Gas grams/mole 

20 n Viscosity of Gas Mpoise 

P Pressure of Gas psia 

p Density of Gas grams/ cm 3 

Q Power of Heat Release Rate watts 

R 0 Resistance at Room Temperature ohms 

25 t Time s 

TABLE II - NOMENCLATURE (con»t) 
Symbol Units 
T Absolute Temperature °C 

U Bridge Output or Amplified Bridge V 

30 Output 

V Volume of Gas or Solid (Microbridge) an 3 

v Speed of Propagation cm/s 

x Temperature coefficient of resistance °c 1 
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SDBSCRIPTS 
c Conduction 
s Microbridge or Solid 

g Gas 

o Room, Reference or Gas Temperature 

Without Microbridge Heating 
h Heater or Hot 

m Middle or Medium 

tl-ti n, A, k or c p at various temperatures 



Figure 4(b) depicts an arrangement wherein the 
center microresistance structure 126 is used as a heater 
flanked by two symmetrically located outer sensing 
resistance elements 122 and 124. The elements 122 and 
15 124 are separated from the heater 126 by a narrow gap. 

Figure 4(c) shows an embodiment configuration in 
which the left element of the bridge 122 is used as the 
heating element and the right element 124 as the sensor. 
This embodiment takes advantage of a rather large 
20 central gap to achieve improved thermal isolation 
between the heater and the sensor. 

Figure 6 shows a modified control circuit which 
uses the center microresistance 126 as heater, while the 
sensing task is performed by the two resistors 122 and 
25 124. The dual heater sensor configuration corresponds 
to Figure 4(b) and the circuit is representative of a 
typical sensor/measurement circuit. Figure 6 includes a 
timer 140 providing square-wave electrical pulses to the 
heater 126. The heater couples the heat pulse to the 
sensors 122 and 124 in the bridge 142. The output of 
the bridge is connected through an amplifier 143 to a 
pair of comparators 144 and 145 which operate "start" 
and "stop" inputs to a counter 146 which counts 10 mHz 
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clock pulses. The counter counts measure the time 
interval (t^ - t t ) between temperatures T 2 and T t 
illustrated in Figure 5; 

Figure 7 is similar to Figure 6 but more detailed. 
5 The bridge configuration is the heater-space-sensor 

configuration of Figure 4(c). The sensor resistance arm 
of the microbridge is set into a wheatstone bridge 150 
at 124. Another proximate resistive arm 122 is fed a 
voltage pulse from pulse generator 151 to provide a heat 

10 pulse into the microbridge element 126. The wheatstone 
bridge 150 also may contain a nulling balancing resistor 
152 which can be used to initially zero the device. The 
microbridge resistor sensor 124 in the wheatstone bridge 
receives the heat pulse from heater element 122 

15 principally by thermal conduction through the 

surrounding fluid. Some conduction, of course, does 
occur through the solid microbridge substrate and 
surroundings. 

The circuitry of Figure 7 is conventional and can 

20 readily be explained with reference to its functional 
operation with regard to processing the bridge output 
signal. The voltage output signals of the bridge 150 
are amplified by differential amplifiers 153 and 154 in 
a differential amplifier section. The imbalance signal 

25 is further amplified by a high gain amplifier at 155. 

The signal at 156 is in the form of a DC voltage signal, 
U, the amplitude of which is solely related to the 
thermal conductivity of the fluid of interest. 

The remainder of the circuitry of Figure 7 includes 

30 a DC level clamping amplifier 157 and isolation 

amplifier 158. The temperature level, time-related 
switching and counting circuitry includes comparators 
159 and 160 together with NAND gates 161 and 162 having 
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outputs which are connected to the counter timing device 
(not shown) as in Figure 6. By measuring the time 
needed for the sensor temperature to rise or fall 
between two or more known temperature values or markers 
5 as represented by sensor resistance or bridge voltage 
outputs, a measure related to the specific heat per unit 
volume, c p , of the fluid of interest is obtained. The 
timing device may be a conventional 10 MHz pulse counter 
or the like. Again, this is illustrated schematically 

10 in Figure 5. 

The output signal from the wheatstone bridge, U, 
represents the voltage imbalance caused by the 
temperature change in microbridge sensor or sensors 
induced by the corresponding heater pulse output. 

15 Because the magnitude of this imbalance is related 

directly to the amount of energy absorbed by the sensor 
or sensors, the amplitude of the signal is directly 
related to the thermal conductivity, k, of the 
conducting media in a manner next explained. 

20 Figure 5 shows that during much of the about 100ms 

wide pulse period, the temperature of the sensor reaches 
and maintains a constant value. During this time, the 
influence of the energy sink or source terms represented 
by specific heat are zero, which means that only thermal 

25 conductivity governs the value of the sensor 
temperature. 

The literature value of the thermal conductivity of 
several gases was compared with the measured sensor 
temperature expressed directly in terms of the measured 
30 wheatstone bridge and balance potential, U. This 

relationship was derived empirically for a microbridge 
of the type depicted in Fig. 4(c) using the least 
squares method in a multiple regression analysis to 
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achieve a best fit curve. This relation can be 
linearized over a modest span sufficient for the purpose 
of the sensor. Other combination configurations of 
heater/ sensor embodiments can likewise be calibrated 
5 using known gases or gases of known k. Thus, using an 
off the shelf flow sensor of the type of Fig. 4(c) in 
the circuit of Fig. 7, a 4volt pulse of 100 ms duration 
was used. 

This yielded an approximate linear relationship 
10 between U and k g of the form: 

k fl = a A U + a 5 (10a) 
The linear approximation holds over enough of a span to 
provide accurate measurements. Similar relations may be 
derived under other measurement conditions including 
15 additional pressure correction terms. 

Further details related to determining the 
coefficients for the algorithms to compute c p are 
described next. This determination requires that the 
measuring system be calibrated first, which consists of 
20 determining the coefficients a 1# a 2 , and a 3 , of the 
algorithm to the computer c p . 

Assuming a two-dimensional model for heat transfer 
in the microbridge, see Figures 4 (a) -4(c), the measured 
sensor temperature response may be described with 
25 reference to the following processes (at zero gas flow) : 

1) Heat release by the heater element film; 

2) Temperature build up in the heater element 
material (FeNi or Pt) and surrounding support 
material (insulator Si 3 NJ , i.e. within the 

30 bridge material; 

3) Conduction towards the sensor via a) the 
bridge material, and b) the fluid phase 
surrounding the bridge; 
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4) Temperature build up in the sensor material 
(as in heater material in item 2 above) , and 
in the gas surrounding it by the heat arriving 
via the above processes; 
5 5) Achieving a steady-state distribution of 

temperature; and 
6) The reverse process to steps 1-5 during the 
start of the heater off -period. 
Further assuming, for the sake of simplicity, 
10 that the specific heats of the involved gaseous and 
solid materials do not depend on temperature, we can 
approximately describe the above processes by the 
following expressions (see Table II above for symbol 
explanation) using the same process numbering as above: 
15 1) Q - V Z /(R 0 (1 + e(T h -T 0 ) ) for small temperature 
rises. 

2) The heater temperature results from balancing the 
heat input and output rates; T h -T 0 - 0/(1^/1^ + 
k s B/L 9 ) with Q in watts; the temperature T h is 

20 established in a time that is short compared to the 

time it takes to reach the sensor if the sensor is 
not identical to the heater, as in configurations 
4 (b) and 4 (c) . 

3) In a truly one-dimensional case most of the 

25 released power Q eventually arrives at the sensor, 

since it only has two ways to go (+x and -x 
directions) . In a two- (or even three-) 
dimensional case a major part of Q gets dissipated 
in the y and z directions, so that only a fraction, 

30 Qp, is conducted to the sensor, with a 

corresponding drop of the original temperature, T h 
down to an intermediate temperature T 0 . The sensor 
then experiences an energy rate arrival of 
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Qc = <V T o> (WL 8 + W L 9> < 1X > 
4) The sensor temperature rise rate is governed by the 

specific heat of the gas surrounding the sensor and 

the closely coupled material of the sensor itself so 

5 that: 

% = (dT/dt) Cps V s + (dT/dt) Cpfl V g . (12) 
It is readily apparent from equation (12) that c„ 
could be determined for an unknown gas if the various 
quantities entering in equations (11) and (12) were 

10 either known or measurable. It has been found, however, 
that even if only dt, dT, T 0 , P and k g are conveniently 
measurable, the other quantities may be determined by 
calibration. This can be done as follows: 

For calibration, gases of known composition 

15 (preferably but not necessarily pure) , and therefore of 
known specific heat and thermal conductivity at the used 
pressure and temperature (both also measured) , are 
brought in contact with the sensor. The effect of the 
pulsed heat releases is recorded in terms of the lapsed 

20 time, t^-t,, as has been described. After noting results 
for various gases, pressures, heater temperatures and/or 
heating/ cooling periods, with pulses of constant 
temperature, voltage, current or power, the recorded 
time and condition data are entered into an array of 

25 data ports which can be used for automatic or 

computerized data processing or other number crunching 
techniques. 

The process can be illustrated with the help of 
equations (11) and (12), by way of example, without 
30 excluding other, similar approaches likely to occur to 
one skilled in numerical analysis. With this in mind, 
the following ports receive data or input for various 
gases, pressures (and temperatures) : 
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Ports: Y XI X2 

Inputs: C M P/P 0 (V*!^ t 2" tl 

A known and available multiple linear regression 
analysis (MLRA, see Figure 8) program can determine the 
5 linear coefficients a,, a^ and a s (e.g., by matrix 

inversion) , which, together with the above input data, 
forms the calibrated expression derived from equations 
(11) and (12) to compute specific heat, c p ; 

c M P/P 0 » a,(Vt,) \ + a 2 (V*i> " a ' (13> 
10 The determined (calibration) coefficients, of 

course, represent the lumped factors of several sensor 
properties or conditions from equations (13) and (14) : 

^ = (T B -TJ(B/L g )(V 8 dT), 

a 2 = (T o -T 0 ) (B g /L s ) (V s dT)k s ( 14 > 

15 ^-^ s /\ 

In order to minimize differences xn T m at the 

sensor location, the most advantageous operation from 

among constant temperature, voltage, current or power is 

chosen. The above method is demonstrated on the basis 

20 of 1) constant voltage pulses, which result in quasi 
square wave heat pulses released by the heater, and 2) 
changes in gas type (CH 4 , Cft, air and 0 2 ) and pressure; 
the chosen configuration was 4 (b) . 

The dt = t^t, and pressure data for each of the 

25 gases used were obtained, for which the c p and k values 
can be obtained from the open literature. This relation 
is linearized by applying the least squares method in a 
multiple linear regression analysis to achieve the best 
fit line. After entering these data into the above 

30 ports Y, XI and X2, the regression analysis program was 
performed. The obtained result was, for a configuration 
as in Figure 4 (b) : 
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a, - -16509, a 2 = 3,5184 and a 3 = .005392 (15) 

The final step in using this calibration method 
involves known means to store, write or burn in the 
5 obtained, tailored values of a 1# a 2 and a 3 for the 

individual microbridge, which may be a Honeywell MICRO- 
SWITCH Model No. AWM-2100V, into the memory linked to 
it. The microsensor is then ready for use to measure 
the specific heat of unknown gases, provided that P and 

10 k be known at the time of measurement. 

Figure 8 depicts a schematic block diagram of a 
device for measuring c p and k. The system includes the 
signal processing circuitry indicated by 170, a jnultiple 
linear regression analysis (MLRA) unit 171 for deriving 

15 the known equation constants for the particular 

microbridge configuration and circuitry used, i.e., a 1 - 
a n , a data bank 172 for storing calibration c p and k 
data and an output interface unit 173. 

With respect to the embodiment of Figure 8, prior 

20 to use, field recalibration may be accomplished simply 
by entering the P, c p and k values of the test gas into 
the data bank. If P cannot be measured independently of 
the sensor already in the subject system its errors can 
be incorporated as a correction in the c p and k 

25 calibration. The measured values of U and dt are then 
used as in the measurement mode to determine sensor 
values of k and c p . If they disagree from the entered 
values the constants a 3 and aj may be modified to fit 
the entered or book values. 

30 This approach may be a practical one for field use, 

but it should be checked by using a second test gas. If 
that agrees, the recalibration may be completed. If 
not, a complete calibration of all a t -a^ coefficients 
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should be made. 

It should be mentioned that in all of the above 
discussion the influence of temperature was not 
mentioned for the sake of simplicity. It is well known, 
5 however, that temperature does influence both c p and k, 
but can be addressed, if necessary, in one of the 
following ways: 

1) Controlled, (expensive and energy consuming) 

or 

10 2) Compensated by special temperature-sensitive 

elements in the analog part of the circuit, or 
3) Entered into the sensor algorithm as an 

additional parameter, which is sensed, 
e.g., by monitoring one of the many available 

15 temperature dependent resistors on the sensor. 

This is the preferred approach for sensing 
systems requiring maximum accuracy. 
With respect to use of the instrument of Figure 8, 
the U and dt = t^t, (and P) signals obtained for an 

20 unknown gas are processed as follows in this mode: 

1) Computation of k from expression (3) using the 
coefficients a 4 and a 5 which have been stored 
in (or burned into) the sensor.' s memory after 
calibration, and 

25 2) Computation of c p from expression (6) . It 

should also be noted that a pressure signal is 
also needed as a basic ingredient since c p is 
used here in relation to a volume of gas as 
opposed to k which is largely pressure 

30 independent if the sensor is used at or above 

atmospheric pressure, at which the gas mean 
free path is small compared to the 
characteristic dimensions of the involved 



WO 93/08457 



PCT/US92/09150 



-37- 

sensor. 

As noted above, the thermal conductivity and 
specific heat microsensor just described comprises 
applicants preferred arrangement for obtaining said 
5 measurements. However, the appended claims should not 
be read so narrowly, but, rather, should be read to 
encompass any known means for determining gas thermal 
conductivity and gas specific heat. 

Now returning to Fig. 3, and as noted above, 

10 several techniques are available in the open literature 
for ascertaining fuel gas viscosity. For example, 
circuitry 88 and sensor 52 could be implemented as a 
flowing gas capillary-type arrangement. (However, 
sensor 52 would need to be moved from area 56 to a 

15 flowing gas environment within chamber 36.) Preferably, 
the novel approach described in the above referenced 
copending application entitled "Multiple Gas Property 
Sensor" is used to ascertain the viscosity of the fuel 
gas. Briefly described, this approach utilizes the 

20 frequency change and series resistance change of a 
crystal resonator positioned within the fuel gas to 
determine both pressure and the combined term molecular 
weight multiplied by viscosity (or, alternatively, 
density multiplied by viscosity) . Using this approach, 

25 sensor circuitry 80 and sensor 46 could be eliminated 
since the process described in said copending 
application produces an absolute pressure reading, which 
is used for correction of measured viscosity and 
molecular weight. For a more complete understanding of 

30 this approach, reference should be made to said case. 

While the invention has been described in detail 
herein in accordance with certain preferred embodiments 
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thereof , many modifications and changes therein may be 
effected by those skilled in the art. Accordingly, it 
is intended by the appended claims to cover all such 
modifications and changes as fall within the true spirit 
5 and scope of the invention. 
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What is claimed is: 

1. A process for the combust ionless measurement of 
the quality of gaseous fuel fed to gas consumption 
devices, said process comprising the steps of: 

(a) conducting at least a partial stream of 
5 the fuel gas through a sensor chamber having a 

first sensor and a second sensor in contact with 
said fuel gas; 

(b) generating a first electrical signal at 
said first sensor, said first electrical signal 

10 being representative of a first fuel gas quality, 

said first fuel gas quality comprising one of the 
thermal conductivity, specific heat, viscosity and 
optical absorption of said fuel gas; 

(c) conducting said first electrical signal to 
15 a computing means; 

(d) generating a second electrical signal at 
said second sensor, said second electrical signal 
being representative of a second fuel gas quality, 
said second fuel gas quality comprising one of said 

20 thermal conductivity, specific heat, viscosity and 

optical absorption fuel gas qualities, said second 
gas quality comprising a different one of said fuel 
gas qualities than said first gas quality; 

(e) conducting said second electrical signal 
25 to the computing means; and 

(f) using said computing means to derive a 
signal of measurement using said received first and 
second electrical signals representative of said 
first fuel gas quality and said second fuel gas 

30 quality, respectively, as a measure for at least 

one of the heat content, density and percent gas 
inerts of the fuel gas according to the 
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formula: 

M = a c + b^ 1 + c,yP 1 

35 where: 

/i « said at least one of said fuel gas heat 

content, density and percent concentration of 
gas inert s, 

a o' b u c v = const^ 5 ' 
40 ml, pi, = exponents, 

x = signal representative of the first fuel gas 

quality, 

y = signal representative of the second fuel gas 
quality, and 

45 + = one of addition, substraction, multiplication 

and division. 

2. The process according to claim 1, wherein said 
fuel gas conducting step (a) is continuous, and further 
comprising the step of automatically periodically 
repeating steps (b)-(f). 

3. The process according to claim 2, further 
comprising the step of: 

converting each derived signal of step (f ) 
into an electrical signal comprising a fuel gas 
5 regulation signal. 

4. The process according to claim 2, further 
comprising the step of: 

transmitting said derived signals of step (f) 
to at least one of display means and recording 
5 means. 

5. The process according to claim 1, wherein said 
first fuel gas quality and said second fuel gas quality 
are sensed at a first temperature, and said process 
further comprises the steps of: 

5 modifying the temperature of at least one of 
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said sensors to a second temperature; 

generating a third electrical signal 
representative of said second fuel gas quality at 
said second temperature; 
10 conducting said third electrical signal to the 

computing means; and 
wherein in said deriving step (f ) , said first, second 
and third electrical signals are used as a measure for 
at least one of the heat content, density and percent 
15 gas inerts of the fuel gas according to the formula: 

where: 

M 555 said fuel gas heat content, density or percent 

concentration of gas inerts, 
a c , bj, Cj, c 2 — constants, 
20 ml, pi, p2 = exponents, 

x t1 = signal representative of the first fuel gas 

quality at temperature tl, 
y t1 = signal representative of the second fuel gas 
quality at temperature tl, and 
25 y t2 = signal representative of the second fuel gas 

quality at temperature t2. 
6. The process according to claim 5, where in said 
signal deriving step (f) formula: 
x = f(n) 

Yti " f i(*ti< k t2> 

5 Ytz - f 2 ( k ti< k t2> 

where: 

n = viscosity of the fuel gas 

= thermal conductivity of the fuel gas at 

said first temperature, tl, 
10 fct2 * thermal conductivity of the fuel gas at 
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said second temperature, t2, 
f (n) = function of fuel gas viscosity , 
f^ktf, k^) « first function of fuel gas thermal 
conductivity at the first 
5 temperature, tl, and thermal 

conductivity at the second 
temperature, t2, and 
f 2 (k t1 , Ic^) = second function of fuel gas thermal 
conductivity at the first 
0 temperature, tl, and thermal 

conductivity at the second 
temperature, t2. 

7. The process according to claim 6, wherein said 
second and third electrical signals are generated by the 
second of said sensors, and said process further 
comprising the step of controllably increasing the 

5 temperature of said second sensor to attain said second 
temperature. 

8. The process according to claim 7, wherein said 
first electrical signal is generated at said first 
temperature. 

9* The process according to claim 8, wherein said 
first temperature is approximately 70°C and said second 
temperature is approximately 120°C. 

10. The process according to claim 6, where in 
said signal deriving step (f) formula: 

f (n) = d/n, d 1 = constant; 
ftCk^k*) - and 

11. The process according to claim 10, wherein the 
fuel gas comprises natural gas and the first, second and 
third electrical signals are used in step (f ) as -a 
measure of the heat content of the gas and wherein: 
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5 a o « 3643.53 

b, * 1050.71 
d, * 102 
ml = 3 

c t * -7.60221 
10 pi - 1 

C 2 » -2294.2 
p2 = 1. 

12. The process according to claim 10 , wherein the 
fuel gas comprises natural gas and the first, second and 
third electrical signals are used in step (f) as a 
measure of the density of the gas and wherein: 

5 a c « 4.3077 

b, « .22937 
d, « 102 
ml = 3 

C t « -.012094 
10 pi = 1 

Cj « -2.2881 
p2 = 1. 

13. The process according to claim 10 , wherein the 
fuel gas comprises natural gas and the first, second and 
third electrical signals are used in step (f) as a 
measure of the percent concentration of inerts including 

5 nitrogen, oxygen and carbon dioxide in the gas and 
wherein: 

a o « 288.69 

b, « -23.818 
d 1 « 102 

10 ml = 3 

c, « -.59575 
pi = 1 

Cj « -173.65 
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p2 - 1. 

14. The process according to claim 10, wherein the 
fuel gas comprises natural gas and the first, second and 
third electrical signals are used in step (f ) as a 
measure of the percent concentration of nitrogen and 

5 oxygen in the gas and wherein: 
a 0 « 464.65 
b t « 9.8185 
d, ~ 102 
ml = 3 

10 c, ~ -.42180 

Pi - 1 

c 2 « -356.18 

p2 = 1. 

15. The process according to claim 10 , wherein the 
fuel gas comprises natural gas and the first, second and 
third electrical signals are used in step (f ) as a 
measure of the percent concentration of carbon dioxide 

5 in the gas and wherein: 
a » 175.96 

O 

b t « -33-636 
dj « 102 
ml = 3 

10 C t « .1739 

pi - 1 

« 182.52 
p2 = 1. 

16. The process according to claim 6, further 
comprising the steps of: 

generating at said at least one sensor a 
fourth electrical signal, said fourth electrical 
5 signal being representative of fuel gas specific 

heat at one of said first and second temperatures ; 
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conducting said fourth electrical signal to 
the computing means; 

using selected ones of said first, second, 
10 third and fourth electrical signals to compute one 

of the molecular weight and density of the fuel 
gas; and 
wherein: 

said step (f) includes deriving said signal 
15 using said first, second, third and fourth 

electrical signals to determine the heat content of 
the fuel gas according to the formula: 
H c - a o + MftnHz)}* 1 + c^Oc^, k t2 )P 1 

i + C 2 f 2< k t1' Kl)* + d 1 C p q1 

20 where: 

H c = heat content of fuel gas; 
a 0 , b v c 1# c 2 , d, = constants; 
ml, pi, p2, ql » exponents; 

z - molecular weight, M, or density, p, of the fuel 
25 gas; 

Cp = specific heat of the fuel gas at one of said 
first and second temper atures. 

17. The process according to claim 16, wherein 
said fuel gas viscosity and said at least one of 
molecular weight and density of said fuel gas are 
determined simultaneously as a combined term viscosity 

5 multiplied by molecular weight of said gas or viscosity 
multiplied by density of said gas. 

18 . The process according to claim 17 , wherein the 
fuel gas comprises natural gas and the first, second, 
third and fourth electrical signals are used in step (f) 
as a measure of the heat content of the gas, and the 

5 combined term viscosity multiplied by molecular weight 
is used in said formula of step (f) and wherein in said 



WO 93/08457 



PCTAJS92/091S0 



-46- 



formula: 



a„ « -1287.7 

O 



b t « -.00090189 
10 ml « 1.7514 

C t * 1,048,800 
pi » -1.7142 

c 2 = 0 

15 p2 = 0 

f z (Kv **) - 0 

dj « 808,700 
ql « .73846. 

19. The process according to claim 1, further 
comprising the step of monitoring and controlling the 
fuel gas temperature within said sensor chamber. 

20. The process according to claim 1, further 
comprising the steps of: 

generating a third electrical signal at one of 
said sensors, said third electrical signal being 
5 representative of a third fuel gas quality, said 

third gas quality comprising a different one of 
said thermal conductivity, specific heat, 
viscosity, and optical absorption fuel gas 
qualities than said first gas quality and said 
10 second gas quality; 

conducting said third electrical signal to 
said computing means; and 
where: 

said signal deriving step (f) uses a formula 
15 of the form: 

11 = ^ + b^ 1 + c^ 1 + djZ* 1 

where: 

fL = said at least one of said fuel gas heat 
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content, density or percent concentration of 
20 gas inert s, 

a 0 , b 1# c 1# d t = constants, 
ml, pi, ql = exponents, 

x « signal representative of said first fuel gas 
quality, 

25 y = signal representative of said second fuel gas 

quality, and 

2 = signal representative of said third fuel gas 
quality. 

21. The process according to claim 1, wherein said 
fuel gas comprises natural gas, and said process further 
comprises the step of obtaining substantially zero gas 
flow within said sensor chamber adjacent said at least 

5 one sensor before said electrical signals are generated. 

22. An arrangement for the combustionless 
measurement of the quality and/or regulation of the 
quantity of gaseous fuel fed to gas consumption devices, 
said arrangement comprising: 

5 a sensor chamber having a plurality of sensors 

therein; 

means for conducting at least a partial stream 
of a fuel gas through the sensor chamber such that 
said gas is in contact with said plurality of. 

10 sensors ; 

first means for generating a first electrical 
signal at one of said plurality of sensors, said 
first electrical signal being representative of a 
first fuel gas quality, said first fuel gas quality 

15 comprising one of thermal conductivity, specific 

heat, viscosity and optical absorption of said fuel 
gas; 

means for transferring said first electrical 
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25 



signal to a computer; 
20 second means for generating a second 

electrical signal at one of said plurality of 
sensors, said second electrical signal being 
representative of a second fuel gas quality, said 
second fuel gas quality comprising one of said 
thermal conductivity, specific heat, viscosity and 
optical absorption fuel gas qualities, said second 
gas quality comprising a different one of said fuel 
gas qualities than said first gas quality; 

means for transferring said second electrical 
30 signal to said computer; and 

means, using said computer, for periodically 
deriving a signal for at least one of measurement 
and regulation utilizing the first and second 
generated electrical signals as a measure for at 
35 least one of the heat content, density and percent 

gas inerts of the fuel gas according to the 
formula: 

M - a 0 + b,x^ + Cl yP 1 

where: 

40 n = said one of said fuel gas heat content, 

density and percent concentration of gas 
inerts, 
a 0 , b t , c t , - constants, 
ml, pi, = exponents, 
45 x = signal representative of the first fuel gas 

quality, 

y = signal representative of the second fuel gas 
quality, and 

+ - one of addition, substraction, multiplication 
50 and division. 
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23. The arrangement according to claim 22, further 
comprising: 

one of a display means and a recording means; 

and 

5 means for transmitting said derived signal to 

said one of display means and recording means. 

24. The arrangement according to claim 22, wherein 
said first fuel gas quality and said second fuel gas 
quality are sensed at a first temperature, and said 
arrangement further comprises: 

5 means for modifying the temperature at at 

least one of said plurality of sensors; 

third means for generating a third electrical 
signal representative of the second fuel gas 
quality at said second temperature; 
10 means for conducting said third electrical 

signal to the computer; and 

said signal deriving means including means for 
using said first, second and third electrical 
signals to determine at least one of the heat 
15 content, density and percent gas inerts of the fuel 

gas according to the formula: 

where: 

/i = said one of said fuel gas heat content, 
density or percent concentration of gas 
20 inerts, 

a c , b t , c 1# c 2 = constants, 
ml, pi, p2 = exponents, 

x t1 = signal representative of the first fuel gas 
quality at temperature tl, 
25 y tl = signal representative of the second fuel gas 
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quality at temperature tl, and 
YtZ = signal representative of the second fuel gas 
quality at temperature t2. 

25. The arrangement according to claim 24, wherein 
in said formula of said signal deriving means: 

5 = f 2 (k^, k t2 ) 

where: 

n = viscosity of the fuel gas, 

= thermal conductivity of the fuel gas at 
said first temperature tl, 
10 ^ = thermal conductivity of the fuel gas at 

said second temperature t2, 
tin) = function of fuel gas viscosity, 
fjCktf* k^) = first function of fuel gas thermal 
conductivity at the first 
15 temperature, tl, and thermal 

conductivity at the second 
temperature, t2, and 
f 2 (k( t , k^) = second function of fuel gas thermal 
conductivity at the first 
20 temperature, tl, and thermal 

conductivity at the second 
temperature, t2. 

26. The arrangement according to claim 25, where 
in said formula of said signal deriving means: 

f (n) « d,/n, d, = constant; 
= k*; and 
5 ^(k«,k2) = k^/kj,. 

27. The arrangement according to claim 26, wherein 
the fuel gas comprises natural gas and the first, second 
and third generated electrical signals are used by said 
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signal deriving means as a measure of the heat content 
5 of the gas and wherein: 
a 0 « 3643.53 
b 1 « 10§0.71 
d, « 102 
ml = 3 

10 C, « -7.60221 

Pi - 1 

« -2294.2 

p2 - 1. 

28. The arrangement according to claim 26 , wherein 
the fuel gas comprises natural gas and the first, second 
and third generated electrical signals are used by said 
signal deriving means as a measure of the density of the 

5 gas and wherein: 

a o « 4.3077 
b, « .22937 
d, « 102 
ml = 3 

10 c 1 « -.02094 

pi - 1 
C 2 a -2.2881 
p2 = 1. 

29. The arrangement according to claim 26, wherein 
the fuel gas comprises natural gas and the first, second 
and third generated electrical signals are used by said 
signal deriving means as a measure of the percent 

5 concentration of inerts including nitrogen, oxygen and 
carbon dioxide in the gas and wherein.: 
a 0 « 288.69 



b 1 « -23.818 
10 ml 



d t « 102 
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C, « -.59575 

pi = 1 

» -173.65 

p2 = 1. 

30. The arrangement according to claim 26, wherein 
the fuel gas comprises natural gas and the first, second 
and third electrical signals are used by said signal 
deriving means as a measure of the percent concentration 

5 of nitrogen and oxygen in the gas and wherein: 
a 0 « 464.65 
b, « 9.8185 

» 102 
ml = 3 

10 C, « -.42180 

pi = 1 
Cg « -356.18 
p2 «= 1. 

31. The arrangement according to claim 26, wherein 
the fuel gas comprises natural gas and the first, second 
and third electrical signals are used by said signal 
deriving means as a measure of the percent concentration 

5 of carbon dioxide in the gas and wherein: 
a ~ -175.96 

O 

b, ~ -33.636 
d, ~ 102 
ml = 3 

10 c, ~ -.1739 

pi = 1 
Cj « 182.52 
p2 = 1. 

32. The arrangement according to claim 22, further 

15 comprising: 

third means for generating a third electrical 
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signal at one of said plurality of sensors, said 
third electrical signal being representative of a 
third fuel gas quality, said third gas quality 

20 comprising a different one of said thermal 

conductivity, specific heat, viscosity and optical 
absorption fuel gas qualities than said first gas 
quality and said second gas quality; 

means for conducting said third electrical 

25 signal to said computing means; and 

where: 

said signal deriving means uses a formula of 
the form: 

fi = a 0 + b^ 1 + c^p 1 + djZ* 1 

30 where: 

/i = said fuel gas heat content, density or percent 

concentration of gas inerts, 
a c , b,, c 1# d, « constants, 
ml, pi, ql = exponents, 
35 x = signal representative of said first fuel gas 

quality, 

y = signal representative of said second fuel gas 
quality, and 

z = signal representative of said third fuel gas 
40 quality. 

33. The arrangement according to claim 22, wherein 
said sensors comprise microsensors and said arrangement 
further comprises means for obtaining substantially zero 
gas flow surrounding said microsensors within said 

5 chamber. 

34. The arrangement according to claim 33, wherein 
said sensor chamber includes a gas inlet and a gas 
outlet and wherein said means for obtaining 
substantially zero gas flow within said chamber adjacent 
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5 said sensors comprises at least one flow restricting 
screen positioned within said sensor chamber between 
said sensor and said gas inlet and outlet. 

35. The arrangement according to claim 34 , wherein 
said first electrical signal generating means and said 
second electrical signal generating means utilize at 
least one mass air flow sensor. 

36. The arrangement according to claim 35, wherein 
said first electrical signal generating means and said 
second electrical signal generating means utilize at 
least one microbridge sensor. 

37. An arrangement for the combustionless 
measurement of the quality and/or regulation of the 
quantity of gaseous fuel fed to gas consumption devices, 
said arrangement comprising: 

5 a sensor chamber having a plurality of sensors 

therein; 

means for conducting at least a partial stream 
of a fuel gas through the sensor chamber such that 
said gas is in contact with said plurality of 
10 sensors; 

first means for generating a first electrical 
signal at one of said plurality of sensors, said 
first electrical signal being representative of one 
of the molecular weight and density of said fuel 
15 gas; 

means for transferring said first electrical 
signal to a computer; 

second means for generating a second 
electrical signal at one of said plurality of 
20 sensors , said second electrical signal being 

representative of the viscosity of said fuel gas; 

means for transferring said second electrical 
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signal to said computer; 

third means for generating a third electrical 

25 signal at one of said plurality of sensors, said 

third electrical signal being representative of the 
thermal conductivity of said fuel gas; 

means for transferring said third electrical 
signal to said computer; 

30 fourth means for generating a fourth 

electrical signal at one of said plurality of 
sensors, said fourth electrical signal being 
representative of the specific heat of said fuel 
gas; 

35 means for transferring said fourth electrical 

signal to said computer; and 

means, using said computer, for periodically 

deriving a signal for at least one of measurement 

and regulation utilizing the first, second, third 
40 and fourth generated electrical signals as a 

measure for at least one of the heat content, 

density and percent gas inerts of the fuel gas 

according to the formula: 

M = a 0 + ^{(nXz)} 01 + c^ 1 + d^ 1 
45 where: 

a Q , b v c 1# d 1 = constants, 

ol, ml, pi = exponents, 

n = viscosity, 

z = molecular weight, M, or density, p, of the fuel 
50 gas, 

k - thermal conductivity, and 
c p = specific heat. 

38. The arrangement according to claim 37, wherein 
the fuel gas comprises natural gas and the first, 
second, third and fourth electrical signals are used by 
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said signal deriving means as a measure of the heat 
5 content of the gas, and the combined term viscosity 

multiplied by fuel gas molecular weight is used in said 
formula, and wherein: 
a « -1287.7 

O 

b t « 00090189 
10 ol * 1.7514 

c t « 1,048,800 
ml « -1.7142 
4, « 808,700 
pi « .73846. 
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